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Abstract. In this talk I review the status of B — B mixing in the Standard Model and the 
room for new physics in B s — B s and Bd — Bd mixing in the light of recent LHCb data. 
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<D • 1. Introduction 



The only source of flavour-changing transition in the Standard Model (SM) is the Yukawa 
interaction of the Higgs doublet H with the fermion fields. The Yukawa lagrangian of the quark 
fields reads 



-L Y = Y d k Q{Hd k R + Yf k Q{Uu k R + h.c. (1) 

00 ' 

Here Q L denotes the SU(2) doublet of the left-handed quark fields of the j-th generation and 
d R , u R are the corresponding right-handed singlet fields. The 3x3 Yukawa matrices Y ,u and 
the Higgs field vacuum expectation value v = 174 GeV combine to the quark mass matrices 
M ,u = Y ' u v. The gauge interactions of the quark fields do not change under unitary rotations 
of any of Q-[, d R , and u R in flavour space. Using these unphysical transformations one can bring 
Y u and Y d to the form 
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Y u = Y u = y c and Y d = V ] Y d with Y d = y s (2) 



and Hi > 0. With the choice adopted in Eq. ([2]) the mass matrix M u of up, charm, and top 
quark is diagonal. The diagonalisation of M requires the additional rotation d ] L = Vjkd^', which 
puts the unitary Cabbibbo-Kobayashi-Maskawa (CKM) matrix V into the W boson vertices: 
WpU^-y^d^ = WfjYjkU^^dj^ . Eq. ([2]) defines a basis of weak eigenstates which is particularly 
suited to display the smallness of flavour violation in the SM. All flavour-changing transitions 
originate from Y d which reads 
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Figure 1. Left: SM box diagram describing B s — B s mixing, a |AB| = 2 process. Right: b — > s 
penguin amplitude, a |Ai?| = 1 transition. 

if all Yukawa couplings are evaluated at the scale n = nit. The off-diagonal element largest 
in magnitude is V t *yb = V^Vb = — 6 • 10 -4 . Moreover, in the SM flavour-changing neutral 
current (FCNC) transitions involve an additional loop suppression, making FCNC processes an 
extremely sensitive probe of physics beyond the SM. This feature puts flavour physics into a 
win-win position: If CMS and ATLAS find new particles, FCNC observables will be used to 
probe the flavour pattern of the BSM theory to which these particles belong. If instead CMS 
and ATLAS do not find any BSM particles, flavour physics will indirectly probe new physics to 
scales exceeding 100 TeV, well beyond the center-of-mass energy of the LHC. 

The weak interaction makes the neutral K ~ sd, D ~ cu, Bd ~ bd and B s ~ bs mesons 
mix with their antiparticles, K, D, Bd, and B s , respectively. This means that, say, a B& meson 
evolves in time into a quantum-mechanical superposition of a B& and a B&. This feature is a 
gold-mine for new-physics searches, as it permits to probe CP phases in any decay to a final 
state which is accessible from both the B& and the B& component of the decaying state. In a 
tree-level decay like B& — > J/tfjKs the time-dependent CP asymmetry probes new physics in the 
mixing amplitude itself. If instead a rare decay (with FCNC decay amplitude) is studied, the 
time-dependent CP asymmetry may reveal a new physics contribution to the decay amplitude. 
Meson-antimeson mixing is a \AF\ = 2 process, meaning that the relevant flavour quantum 
numbers change by two units. In the case of B s — B s mixing depicted in Fig. [TJ these flavour 
quantum numbers are beauty B and strangeness S. Weak decays are |A.F| = 1 processes, 
Fig. [1] shows a \AB\ = \AS\ = 1 FCNC decay amplitude. It is instructive to compare the 
reach of \AF\ = 2 and \AF\ = 1 amplitudes to new physics: The flavour violation in the SM 
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the amplitudes ^4g M and A SM in Fig. [I] is governed by the small CKM element Vt s and 
both diagrams scale as 1/M^. A new-physics (NP) contribution may instead involve some new 
flavour-violating parameter <5fcnc- If ^ ne sca le of NP is A and the NP contribution enters at 
the one-loop level, one finds 

, ,|A_B|=2, , j. |0 ,,9 I a\ ab \ =1 \ \s i t\,t2 
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One sees that B s — B s mixing is more sensitive to generic NP than an FCNC b — > s decay 
amplitude, if |<5fcnc| > Wts\ ~ 0.04. The estimate in Eq. (JH) applies, for example, to the 
(N)MSSM with generic flavour structure in the bilinear SUSY-breaking terms, which induce new 
bfi — >■ sr or bz — > sl transitions through squark-gluino loops [HE]. But \AB\ = 1 transitions 
can be more sensitive probes of NP in models in which 6r — > sl or bi — > sr transitions are 
parametrically enhanced over their SM counterparts. This situation occurs in the (N)MSSM 
with a large value of the parameter tan/3 or large trilinear SUSY-breaking terms. Theories 
with such chirally enhanced FCNC transitions are efficiently probed with the radiative decays 
b — > sj, dj [3ttT0]. leptonic decays like the recently observed [UJ decay B s — >• u^ pT [71 [±21421] 
and the chromomagnetic contribution to non-leptonic decays like Bd — > 4>Ks [18J. 



2. B — B mixing: formalism and Standard-Model prediction 

B q — B q mixing with q = d or q = s is governed by the 2x2 matrix M q — iT q /2, with the 
hermitian mass and decay matrices M q and T q . Due to non- vanishing off-diagonal elements 
M q 2 and Y\ 2 a meson tagged as a B q at time t = will for t > evolve into a superposition of 
B q and i? 9 . In the SM M[ 2 is calculated from the dispersive part of the box diagram in Fig. [H 
which amounts to discarding the imaginary part of the loop integral. Tf 2 is instead obtained 
from the absorptive part, meaning that the real part of the loop integral is dropped. Mf 2 and 
r^ 2 are calculated in the same way from the box diagram with the external s quarks replaced by 
d quarks. M q 2 is dominated by the top contribution, while Y\ 2 stems from the diagrams with 
internal up and charm. T\ 2 is made up of all final states into which both B q and B q can decay. 
There are three physical quantities in B q — B q mixing: 



|M? 2 |, |r? 2 |, ^arg(-^p) 



By diagonalising M—i T/2 one finds the two mass eigenstates B q and B^ , where the superscripts 
stand for "heavy" and "light" . These eigenstates differ in their mass and width with 



Am, = M q H -M q L ~ 2|Mf 2 |, Ar g = T q L - T q H a 2|F 
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The mass difference Am q simply equals the frequency at which B q and B q oscillate into each 
other. The width difference Ar s is sizable, so that in general untagged B s decays are governed 
by the sum of two exponentials. There is no useful data on the tiny width difference AI^ in the 
Bd system yet, which is predicted to be around 0.5% of T d L H [2"2Tf2"5] . The CP-violating phase 
4> q can be determined by measuring the CP asymmetry in flavour-specific decays: 
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A decay B q — > f is called flavour-specific, if the decay B q — >• / is forbidden. The standard 
method to determine af s uses semileptonic decays, so that a q s is often called semileptonic CP 
asymmetry. While this measurement simply requires to count the numbers of positive and 
negative leptons in B q decays, it is still very difficult, because |r^ 2 | <^. \M q 2 \ renders |a£ s | small, 
even if NP contributions to M q 2 enhance 4> q over its small SM value. 

2.1. Mass difference Am s 

The theoretical prediction of Am q requires the separation of short-distance and long-distance 



QCD effects. To this end one employs an operator product expansion, which factorises M? 2 as 



M q 12 = (V tq V tb f C(B q \Q\B q ). (5) 

The Wilson coefficient C comprises the short-distance physics, with all dependence on the heavy 
particle masses. QCD corrections to C have been calculated reliably in perturbation theory |26j. 
Since the CKM elements are factored out in Eq. i^, C = C(mt,M\v,ot s ) is real in the SM. 

Q = Q L 'y l/ b L q L j' / b L (6) 

is a local four-quark operator describing a point-like \AB\ = 2 transition, see Fig. [2j The 
formalism is the same for the B^ and B s complex; for definiteness I only discuss the case q = s 
in the following. The average of the CDF and LHCb measurements of Am s ~ 2|Mf 2 | is |27j 

Amf p = (17.719 ± 0.043) ps" 1 , (7) 




Figure 2. The |AJ5| = 2 operator Q of Eq. Q. It is pictorially obtained by shrinking the 
B q — B q mixing box diagram to a point. 



meaning that \Mf 2 \ is known very precisely. However, we need the hadronic matrix element 
(B S \Q\B S ) to confront Eq. © with the SM. The latter is usually parametrised as 

(B S \Q\B S ) = 2 -M 2 B J 2 Bs B Bs . 

Here Mb s and fs s are mass and decay constant of the B s meson, respectively. The hadronic 
parameter Bb s depends on the renormalisation scheme and scale used to define Q, in this talk 
Bb s ~ 0.85 is understood to be evaluated at the scale /i = mj, in the MS scheme. Noting that 
| Vi S | is fixed by CKM unitarity from the well-measured element \V c b\, we can write 



Am s = 18.8 ± 0.6y cb ± 0.3^ ± 0.1 J ps" 1 / ^ B n \ l[ ^ 2 , (8) 



(220 MeV) 5 



where the uncertainties from the experimental errors of the input parameters indicated. 
Averaging various calculations from lattice gauge theory one finds 



f B B Bs = [(211 ±9)MeV] 2 . (9) 

Inserting this result into Eq. (JSJ) implies 

Am s = (17.3 ±1.5) ps" 1 (10) 

complying excellently with the experimental result in Eq. ([7]). However, the average in Eq. ([9]) 
mainly involves different calculations of fs s , which are combined with two fairly old results for 
Bb s - With the recent preliminary Fermilab/MILC result [29j 



one finds instead 



/| b Bb = 0.0559(68) GeV 2 ~ [(237 ± 14) MeV] 1 



Am s = (21.7 ±2.6) ps" 1 . (11) 



Therefore more effort on lattice-QCD calculations of f B B Bs is highly desirable. As long 
as modern calculations of this quantity are unavailable, I recommend to inflate the error in 
Eq. ()10p to 2.5 ps™ 1 . One concludes that the precise measurement in Eq. ([7|) still permits a NP 
contribution of 15% in Am s . 

2.2. CP phase from B s -> J/ ijxj) 

While Am s fixes the magnitude of Mf 2 , the phase of Mf 2 can be probed through the mixing- 
induced CP asymmetry in B s — > J/ipcj). The final state contains two vector mesons, so that the 
orbital angular momentum L can assume the values 0,1, and 2. The final states (J/4' ( I ) )l=q,2 




T—J/ip 



c 



Figure 3. Amplitudes of B s — > J/tp <\> and B s — >• J/ip <p>. The cross denotes the T^-mediated 
b — > ccs decay. 



have CP quantum number rjcp = 1, while (J/iP4>)l=i is CP-odd. Denoting a meson which 
was a B s at time t = by B s (t) (with an analogous definition of B s (t)), one can define the 
time-dependent CP asymmetry 



acpit) 



T(B s (t) -» (J/j;<l>) L )-T(B a (t) -» (J/j^U) 
r(£ s (t) -»• ( J/V 0)l) + r(S,(t) -+ (J/V </>)l) ' 



(12) 



The two interfering amplitudes giving rise to acpit) are shown in Fig. The analytical result 
reads 



acp(t) = rjcp 



sin(2/3 s ) sin(Am s £) 



cosh(Ar s t/2) - cos(2/3 s ) sinh(Ar s t/2) ' 



(13) 



The CP phase entering acp{t) is j3 s = arg(— V t * b Vt s / (V* b V cs )) , in the standard phase convention 
of the CKM matrix (3 S is essentially just the phase of — Vj s . It should be mentioned that Eq. (|12p 
is a theoretical definition, in practice the experimental separation of the L = 0, 1,2 amplitudes 
involves a complicated angular analysis of the J/tp and decays. 

The CKM matrix can be parametrised in terms of the four Wolfenstein parameters \,A,~p, and 
rj |30p31| . Expanding to leading order in A = 0.225 one finds sin(2/3 s ) ~ 2\ 2 rj. The parameter r\ is 
the height of the CKM unitarity triangle defined by ~p + irj = — V* b V u d/ (V* b V c d) ■ The suppression 
by A 2 renders sin(2/3 s ) small; determining the CKM elements from a global fit to the data gives 
rj = 0.343 ± 0.015 and leads to the SM prediction [28] 



2/3 s 



2.1° ±0.1°. 



(14) 



The combination of CDF, DO, ATLAS, and LHCb data on B s — >■ J/ip<p and of LHCb data on 
B s ->• J/ipTT+ir- gives [27] 



2^fP = 0.7 o +4:s°- 



(15) 



Of course, this average is dominated by the LHCb data. At this conference Jeroen van Leerdam 
has reported 2/3 s = 0.1° ±4.8° tat ± 1.5° yst from the combined LHCb analysis of B s ->• J/ipK + K~ 
and B s —> J/ipir + ir~ |32j. In Eq. (|13|) an additional small contribution, the "penguin pollution" 
has been neglected. This -presently uncalculable- contribution inflicts an additional error of 
order 1° on 2/3f xp in Eq. JEJ. 



2.3. Decay matrix T\ 2 

The calculation of T^ 2 , q = d,s, is needed for the width difference Ar g ~ 2|r^ 2 | cos 4> q and the 

semileptonic CP asymmetry aj? s = Ljq , sin^ g . The theoretical results of Refs. [22,231(33] have 

recently been updated to 



</> a = 0.22° ±0.06°, d = -4.3°±1.4°, 

af s = (1.9 ±0.3) -10" 5 , and af s = -(4.1 ± 0.6) • 10 -4 . (16) 

The prediction of T\ 2 involves new operators in addition to Q in Eq. ([6]). However, the matrix 
element of Q comes with the largest coefficient, so that the ratio Ar g / Am, = |r| 2 |/[Af' 2 | suffers 
from smaller hadronic uncertainties than Ar g . Predicting Ar q with the help of the experimental 
values Am^ p = 0.507 ps _1 and Ams Xp in Eq. (|7|) one infers from Ref. [25j : 

AT 

AT d = ^Amf p = (27±5)-W~\s- 1 

AT 
Ar s = — ^Amf p = (0.090 ±0.018) ps" 1 (17) 

Am s 

The more recent update in Ref. |34| differs from Ref. |25j in two respects: The quoted results 
are obtained in a different renormalisation scheme and use the preliminary lattice results of [29] 
instead of the averages from [28] . Ref. [34] finds for Ar s : 

Af 

Ar s = — —Amf p = (0.078 ± 0.018) ps" 1 , (18) 

Am s 

which is consistent with Eq. (|17p . In Eqs. (|17p and (|18p all uncertainties are added in quadrature, 
which may not be a conservative estimate of the overall error. The ranges comply with the LHCb 
measurement 132 



Ar LHCb = [an6 ± .018 stat ± 0.006 syst ] ps- 1 (19) 

and the world average [27] : 

Arf p = [0.089 ± 0.012] ps" 1 . (20) 

3. New physics in B — B mixing 

3.1. A model-independent analysis 

If some NP amplitude adds to the B q — B q mixing box diagram, both \Mf 2 \ and 4> q may deviate 

from their SM predictions. The D0 experiment has measured 



A%£ = (0.532 ± 0.039)4 ± (0-468 ± 0.039)4 

= (-7.87 ± 1.72 ± 0.93) • 10~ 3 , (21) 

which is 3.9<7 off the SM prediction inferred from Eq. (|16|) . 

A SL = (-0.20 ± 0.03) -10" 3 . (22) 

Adding a NP contribution (fi ds to either <^ or <p s in Eq. f)16|) may enhance |^4gL|, but 
the very same contribution will also affect the value of (3 S in Eq. (|15p an d the value of 
(3 = arg = &rg(—VtbV* d / (V c bV* d )) found from acp{t)(B ( i — > J/ipKs)- One can parametrise 
NP in Bd — Bd mixing and B s — B s mixing by two complex parameters A d and A s : 

M q = M SM ' q ■ A A = IA Ip^? 
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Figure 4. Allowed range for A s , see 



for details. 



In summer 2010, before the advent of precision data from LHCb, a global fit of all relevant 
flavour data to the CKM elements and A^ and A s has resulted in a 3.6a evidence of NP, with 
a large negative NP phase </>^ [38]. In spring 2012 this fit has been repeated [39] : Figs. 0] and 
[5] show an update of the results in Ref. [39J with the data of 2012 summer conferences [28] . 
The fit tries to accomodate A^l in Eq. (|2ip and a slightly high value of the world average for 
B(B ->• rv) through </ft < 0. Since a C p(t)(B d ->• J/if)K s ) precisely fixes 2/3 + ^ = 42.8° ±1.6°, 
<pj < comes with a higher value of /3 than in the SM. (/3 > /3 entails a larger value of 
\V u b\ which governs B{B — >• tu).) Contrary to the situation in 2010, the Standard Model point 
A s = Ad = 1 is merely disfavoured by 1 standard deviation, consistent with natural statistical 
fluctuations. At the best-fit point the problem with ^4gL is only marginally alleviated. 

One could relax the problem with ^4gL without affecting acp{t){Bd — >• J/tpKs) and 
a Cp(*)(-Bs - ^ J/^xp) by postulating new physics in Tf 2 or Tf 2 [391 - I4T] . Since Tf 2 originates 
from Cabibbo-favoured tree-level decays, it can hardly be changed in a significant way without 
spoiling the value of the average decay width T s = (T s H +T s L )/2 [551I3T?] . The LHCb measurement 
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Figure 5. Allowed range for A^, see 



for details. 



r LHCb = (0.6580 ± 0.0054 ± 0.0066) ps" 1 [32] implies 



r. 



T B d 



0.997 ±0.013 



(23) 



in excellent agreement with the SM prediction T B Jr Bd = 0.998 ±0.003 (25HI2] . NP in the doubly 
Cabibbo-suppressed quantity Tf 2 is phenomenologically only poorly constrained, but requires a 
somewhat contrived model of NP. 



3.2. A supersymmetric SO(10) model 

The Minimal Supersymmetric Standard Model (MSSM) has many new sources of flavour 
violation, which all reside in the super symmetry-breaking sector. It is no problem to get 
a big effect in a chosen FCNC process, but rather to suppress big effects elsewhere. This 
supersymmetric flavour problem is substantially alleviated with the lower bounds on the squark 
massed placed by ATLAS and CMS. Grand Unified Theories (GUT) offer the possibility to 



have "controlled" deviations from the CKM pattern of flavour violation in the quark sector: In 
GUTs quarks and leptons are unified in common symmetry multiplets, opening the possibility 
to observe the large lepton-flavour mixing encoded in the Pontecorvo-Maki-Nakagawa-Sakata 



(PMNS) matrix C/pmns i n quark flavour physics 



Consider SU(5) multiplets: 
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If the observed large atmospheric neutrino mixing angle stems from a rotation of 52 and 
53 in flavour space, it will induce a large 6r — s/j-mixing. Contrary to the situation with 
right-handed quark fields, rotations of right-handed squark fields in flavour space are physical 
because of the soft SUSY-breaking terms. The Chang-Masiero-Murayama (CMM) model has 
implemented this idea in a GUT based on the symmetry breaking chain SO(10)— » SU(5) — > 



SU(3)xSU(2) L > 
diagonalised as 



U(l)y |45H47| . In some weak basis the Yukawa matrix of down (s)quarks is 



Y d = V ( 



CKM 




C^PMNS 



and the right-handed down squark mass matrix has the following diagonal form: 



m 2 r(M z ) = di&g 



m % 



m 7 , m 7 — At 

cr d a 



The real parameter Aj is calculated from renormalisation-group effects driven by the top- Yukawa 
coupling. 



Rotating Y^ to diagonal form puts the large atmospheric neutrino mixing angle into m^: 
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(24) 



•B s mixing! 



The CP phase £ affects CP violation in B s 

In Ref. [Uj we have confronted the CMM model with flavour data. The analysis involves 
seven parameters, which we have fixed by choosing values for the squark masses Mq, Mr of 
right-handed up and down squarks, the trilinear term af of the first generation, the gluino 
mass 77ip 3 , tan/3, and the sign of the higgsino mass parameter /i. We have considered B s — B s 
mixing, b — > 57, r — > /ry, vacuum stability bounds, lower bounds on sparticle masses and 
the lower bound on the lightest Higgs boson. From these inputs first universal SUSY-breaking 
terms defined at a fundamental scale near the Planck scale have been determined through the 
renormalisation group equations (RGE). Subsequently we have used the RGE to determine all 
low-energy parameters, with the MSSM as the low-energy theory. 

Two experimental results of the year 2012 put the CMM model under pressure: First, the 
sizable neutrino mixing angle #13 leads to an unduly large effect in B(fi — > ej). In Eq. (|24p I have 
tacitly assumed a C/pmns with tri-bimaximal mixing, corresponding to #13 = 0. For the actual 
value 0i3 ~ 8° the (1,2) element of the charged slepton mass matrix gets large, and one has to 
resort to much larger sfermion masses than those considered in Ref. |47j . Second, the Higgs mass 
of 126 GeV challenges the model and we are unable to find parameters which simultaneously 
satisfy the Higgs mass constraint and the experimental upper bound on B(fi — >■ cy) [48J. The 
problem with the Higgs mass could be circumvented by considering the NMSSM as the low- 
energy theory. 



4. Conclusions 

LHCb has provided us with a significantly better insight into the B s — B s mixing complex. 
Am s and AT S comply with the SM, but we need better lattice data for the hadronic matrix 
elements involved. Theoretical uncertainties still permit an 0(20%) NP contribution to the 
B s — B s amplitude M± 2 - While in 2010 the D0 result for Asl could be explained in scenarios 
with NP only in M 12 , the LHCb data on B s — > J /ijxj) now prohibit this solution. An alternative 
explanation invoking new physics in r| 2 is not viable, because this will spoil the ratio T s /T(i 
which agrees well with the SM prediction. Maybe it is worthwile to look at NP in Tf 2 , although 
this possibility leads to somewhat contrived models. Models of GUT flavour physics with b — > s 
transition driven by the atmospheric neutrino mixing angle are under pressure from the large 
value of #i3, which induces a too large B(fj, — > cy). Moreover, in the studied CMM model it 
seems impossible to accomodate the measured value of the lightest Higgs mass, if one insists on 
the MSSM as the low-energy theory. 
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